Zero-bias anomaly in the tunneling density of states of graphene 
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In the vicinity of the Fermi energy, the band structure of graphene is well described by a Dirac 
equation. Impurities will generally induce both a scalar potential as well as a (fictitious) gauge field 
acting on the Dirac fermions. We show that the angular dependence of the zero-bias anomaly in the 
spatially resolved tunneling density of states (TDOS) around a particular impurity allows one to 
distinguish between these two contributions. Our predictions can be tested in scanning-tunneling- 
microscopy measurements on graphene. 
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Introduction — Since its recent experimental realization 
P, Q , graphene - a monolayer of graphite - has attracted 
a lot of attention due to its remarkable electronic struc- 
ture. At low doping, the Fermi surface of graphene lies 
in the vicinity of two points in the Brillouin zone (termed 
Dirac points or valleys) near which the spectrum is char- 
acterized by a massless Dirac dispersion 0, 0] with veloc- 
ity v. Pioneering experiments on this novel 2D electron 
system have shown that the Dirac nature of carriers in- 
duces an anomalous integer quantum Hall effect as well 
as a finite conductivity at vanishing carrier density 0, Q . 

Within the independent-electron approximation, the 
density of carrier states v(E) = \E\/27rh 2 v 2 vanishes lin- 
early at the Dirac point E = 0. Thus, Coulomb inter- 
actions remain unscreened when the Fermi energy is at 
the Dirac point, which leads to deviations from conven- 
tional Fermi- liquid expectations The Fock diagram 
for the electronic self energy, shown in Fig. 1(a), entails 
a logarithmic correction to the linear dispersion relation, 
E(p) = [v + (e 2 /4)ln(D/vp)]p, with D the bandwidth, 
and hence to the tunneling density of states 
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At the same time, there is no renormalization of the elec- 
tron charge from the diagram in Fig. 1(b). Several au- 
thors have discussed possible instabilities of the electron 
system in graphene 0, 0? 01? although at present, there 
is no related experimental evidence. 

Impurities and defects generally induce both a scalar 
potential as well as a fictitious gauge field in the ef- 
fective Dirac description [9]. Site disorder in the un- 
derlying tight-binding model is associated with a ran- 
dom scalar potential. Random hopping due to lattice 
deformations leads to a random gauge field which is 
abelian (nonabelian) for intravalley (intervalley) scat- 
tering. The presence of a random gauge field leads to 
rich weak-localization physics [Io|, EH E^- Several re- 
cent works have addressed localization for Dirac electrons 
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Unlike conventional two-dimensional electron systems 
(2DES), the 2DES in graphene is exposed at the sur- 
face, making it directly accessible to local-probe measure- 
ments, such as scanning tunneling microscopy (STM) or 
scanning single-electron transistors [16[. Motivated by 
this fact, we investigate the tunneling density of states 
(TDOS) in this paper. Unlike the logarithmic corrections 
due to Coulomb interactions which become singular at 
the Dirac point, the combined effects of disorder and in- 
teractions at finite doping lead to a zero-bias anomaly 
(ZBA) which is tied to the Fermi energy. 

The zero-bias anomaly in disordered conductors arises 
from scattering of electrons on impurities and on the po- 
tential generated around the impurity by the Friedel os- 
cillations 
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181 ]. This combined effect of disorder and 
interactions leads to a suppression of the tunneling den- 
sity of states at low biases u. The ZBA in graphene 
raises several issues, (i) The wavelength of Friedel oscil- 
lations in graphene [19[ diverges as the Fermi energy ap- 
proaches the Dirac point. At the same time, the strength 
of the Coulomb interaction increases due to the absence 
of screening at the Dirac point. We study how this com- 
petition affects the magnitude of the ZBA as the Fermi 
energy Ep approaches the Dirac point within the regime 
E F r > 1 (with r the elastic mean free time due to impu- 
rity scattering), (ii) We show that the zero-bias anomaly 
provides a means of distinguishing experimentally be- 
tween the scalar potential and the fictitious gauge po- 
tential induced in the Dirac equation by a specific im- 
purity. Our considerations focus on the quasiballistic 
regime Ep > oo > 1/r relevant in relatively clean sam- 
ples. In this regime, the physics is essentially captured 
by considering electrons scattering off isolated impurities. 
Our approach breaks down at the lowest biases uo < 1/r, 
where electron diffusion becomes relevant [2o| . 

The use of the Dirac formalism restricts us to a dis- 
cussion of the local TDOS which is coarse-grained over a 
region large compared to the lattice spacing. At the same 
time, the long- wavelength nature of the Dirac description 
allows us to include the effects of electron-electron inter- 
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FIG. 1: (a) Fock self energy of Dirac fermions; (b) leading 
vertex correction; (c) Hartree and (d) Fock diagram for the 
zero-bias anomaly. 



actions on the local TDOS which have been neglected 
in several earlier studies performed in the framework 



of a tight-binding model [21|, [22[. As we argue below, 



the effects of electron-electron interactions on the coarse 
grained TDOS are of particular importance in graphene 
since the effects of impurity scattering alone are sup- 
pressed by the chiral symmetry of Dirac fermions. 

Graphene — Within the tight-binding approximation, 
the bandstructure of graphene is described by the Hamil- 
tonian H = —t^2^\ c\cj + h.c. on a hexagonal lattice, 
where t is the hopping matrix element, q is the anni- 
hilation operator of an electron on lattice site z, and 
only nearest-neighbor hopping has been considered. The 
2D hexagonal lattice consists of two identical sublat- 
tices A and B, and thus two sites per unit cell. We 
choose the vectors connecting a B site with the neigh- 
boring A sites as e 1 = a (—1,0), e 2 = a (l/2,\/3/2) 
and e 3 = a (1/2, —y/3/2) with a being the bond length. 
The Bloch spectrum of the tight-binding Hamiltonian 
has zero energy (corresponding to the Fermi energy at 
half filling) at two inequivalent Dirac points in the re- 
ciprocal lattice, which we choose to be k ± = ±k 
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with k D = 27r/(3\/3a) (\/3, 1). Linearizing the spectrum 
about both Dirac points and arranging the Hamiltonian 
in a 4 x 4 matrix form, one can write (see e.g. Ref. [Hj]) 



H = hv £ • k 



(2) 



with v = 3ta/2, k the 2D wavenumber deviation from 
the Dirac point, ^ XjV = H Z <8> cr x y the components of the 
vector T, where 11^ and <jj are Pauli matrices acting in 
the spaces of the Dirac points (+/— ) and the sublattices 
(A/B), respectively. The Hamiltonian in Eq. ([2j) acts 
on the four-component spinors (u\ k? k> U B k' U A k) °f 
Bloch amplitudes in the A/B and +/— channels. 

Within the Dirac Hamiltonian Eq. ([2]), an impurity 
localized, say, at the origin can be quite generally ac- 
counted for by adding a general local operator consistent 
with time-reversal symmetry, 



H dis (r)=uS(r). 
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Here, u is a matrix in the 4x4 Dirac space, reflect- 
ing inhomogeneities in both site energies and hopping 
amplitudes of the underlying tight-binding Hamiltonian 
The matrix u can be parametrized in terms of 



ten real numbers u t 



0? 



(s, I = x,y, z) as 



s,l=x,y,z 



'sl^s A l 



with I = ll n 



(4) 
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^x,y ^~x,y 



cr the identity matrix and Y> z 

The parameters uo,u s i can be classified according to 
whether they describe intravalley or intervalley scattering 
and whether they take the form of potential or (fictitious) 
gauge field disorder in the Dirac equation. For example, 
intervalley scattering will involve H x ,y and contributions 
to the fictitious gauge field involve ^ Xy y The complete 
classification is detailed in Table I. 





potential 


(fictitious) gauge field 


intra-valley 


Uq , U zz 


u xz 5 u yz (abelian) 


inter-valley 


^zx i ^zy 


u xx , u xy , u yx , u yy ( iioiiab elian) 



TABLE I: Physical significance of the various components 
of the disorder potential u in terms of type of potential 
(columns) and type of scattering (rows). 

For Dirac fermions described by Eq. ([2j), the re- 
tarded Green function with energy e and momentum 
p is Gf(p) = ^±s Tp /[e±vp + irj] where s p = 
1/2 [1 + 5] • p/p] denotes the projector onto the state 
with positive chirality XI • p/p and n — > + . In view 
of particle-hole symmetry, we specify attention to elec- 
tron doping (positive Fermi energy). Then, in the regime 
k F r ^> 1, the real space Green function takes the form 



G?(r,0) 



4p e r r 
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with p e = e/Tiv. Here, we retained the next-to-leading 
term in in the square bracket for later convenience. 

The projector s r is defined in analogy with s p . 

Tunneling density of states — To compute the (coarse- 
grained) local TDOS 
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in the quasiballistic regime, we consider contributions to 
the full single-particle Green function due to interac- 
tions and scattering off isolated impurities. The leading- 
order correction arises from paths involving a single im- 
purity scattering event, corresponding to direct backscat- 
tering of electrons to the tunneling point. 

In conventional 2DES, this contribution falls off as 1/r 
and oscillates with wavevector 2kp (r denotes the dis- 
tance of the tunneling point from the impurity). For 
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Dirac electrons, the leading term is suppressed by chi- 
rality and a straight-forward calculation yields a faster 
spatial decay, 



5v F (r) 



-4z/ F u 



sin(2/c F r) 
(k F r) 2 
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Due to this suppression, one expects the combined ef- 
fects of disorder and interactions, involving scattering on 
the impurity as well as on the potential generated by the 
surrounding Friedel oscillations, to be particularly im- 
portant in graphene. It is the latter contribution which 
we now address. 

In the presence of a finite Fermi surface, the impurity 
potential generates Friedel oscillations in the carrier den- 
sity. This in turn yields an additional oscillatory scatter- 
ing potential due to the Coulomb interaction, affecting 
the electronic return probability and consequently the 
TDOS. To first order in the electron-electron interaction 
and the impurity potential, the relevant correction to the 
Dirac fermion Green function [18[ is given by 

50?(r,r) = J d ri dr 2 Gf (r, 0)fiGf (0, r x ) 

xH HF (r l5 r 2 )Gf(r 2 ,r) + (u <- iT HF ), (8) 

in terms of the Hartree and Fock potentials 
#HF(ri,r 2 ) = Vjy(ri)J(ri - r 2 ) - V F (r 1 ,r 2 ) 
with VJy(n) = Tr[/drT(ri - r')fy(r',r')] and 

^F( r l^ r 2) = \ V ( Y 1 ~ r 2)M r l> r 2)' Here > 

p(r, r') = 2Y, k <k F ^k( r ) ® ^k( r ') denotes the den- 
sity matrix and V(r) the screened Coulomb interaction. 
The shorthand (u <-> Hup) denotes a similar contri- 
bution with the order of u and H^f interchanged and 
appropriate changes to the spatial arguments. To first 
order in the impurity potential Eq. (|3|), the density 
matrix emerges from the corrections 
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FIG. 2: Schematic of the stationary paths for the Fock con- 
tribution 



the Hartree contribution to the zero-bias anomaly is sup- 
pressed and to leading order, we only need to consider 
the Fock contribution. Then, the contribution Eq. J5J) 
corresponds to a particle being injected at r, moving to 
the impurity (located at the origin), being scattered to 
r 1 where it experiences the non-local Fock potential, and 
finally returning from r 2 to the injection point r. The 
dominant contribution to the integral can be identified 
by analyzing the phase factors of the integrand, cf. Fig. [2] 
and Ref. [l8[ . Clearly, paths entailing phases which oscil- 
late rapidly with r 1 and r 2 lead to negligible corrections. 
The three retarded Green functions in Eq. J5J) yield a 
phase factor ex.p[ip e (r + r x + |r 2 — r|)] while the density 
matrix yields ex.p[iak F (r 1 + /3r 2 )], with a, /3 = ±1. We 
are interested in corrections to the TDOS in the vicinity 
of the Fermi energy, so that e = E F + fiw, with Tiuj <C E F 



and p e 



uj/v. A correction to the TDOS which 



dp 



varies slowly with r arises from the region r 2 ^> r, so 
that |r 2 — r| ~ r 2 — rcos#, where is the angle be- 
tween r and r 2 . Indeed, the oscillatory dependence on r 
is essentially cancelled if \6\ < 1/ \/p e r ~ l/^/k F r <C 1, 
yielding r 2 ~ r 2 r in leading order. The remaining phase 
dependence on r x and r 2 is exp[zp e (r 1 + r 2 )], which 
can be cancelled by choosing a = — 1 and f3 = 1 in 
the density matrix, i.e., by the rapidly oscillating term 
exp[— ik F (r 1 +r 2 )] in it^(r 2 , r x ). This cancellation is valid 



(2tt) E k -vp + irj 



e ip r u |k) . (9) as long as r 2 < v/cj. As a result, we obtain 



to the plane- wave wavefunctions ^k( r ) = e * k r |k) with 
energy E k . The corresponding correction to the density 
matrix takes the form 5p(r, r') = i?(r, r') -\-R^(r\ r) with 

ik e ik F (r-r') 
e ik F (r+r') 

+ — 4s r us_ r ,] (10) 

in leading order in l/k F r <C 1. 

It has been shown that due to the chiral symmetry 
of Dirac fermions, the Friedel oscillations of the electron 
density decay as 1/r 3 as opposed to 1/r 2 in conventional 
two-dimensional electron systems [19[. For this reason, 



<S „ US „ Uj S y. 



(11) 



with V(0) being the Fourier transform of the screened 
Coulomb interaction at zero momentum. The correc- 
tion to the local TDOS, normalized to the bare DOS 
of graphene at the Fermi level, v F = k F /(27rhv), is then 



<Si/ w (r) ^ 4*4^(0) 



r 3/2 



— Tr [s r us_ 

(k F r) 2 



(12) 



valid in the regime l/k F < r < v/uj. The correction 
saturates for r < l/k F and decays as 1/r 3 , for r > v/uo. 
For a finite density rii of impurities, spatial averaging of 
Eq. ([T2]) (or, alternatively, the diagrams in Figs. 1(c) and 



(d)) yields the average TDOS 
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8v 2 



Tr 



2u 2 



a—x,y 
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where we employed the Thomas- Fermi expression V'(O) = 
\jv F for the 2D Coulomb interaction. Similar to conven- 
tional two-dimensional electron systems [Hj, graphene 
also exhibits a logarithmic zero-bias anomaly. Eq. (p~3|) 
shows that the relative strength of the ZBA is only 
logarithmically dependent on the Fermi energy (since 
v~f oc kp)t even though the wavelength of the Friedel 
oscillations diverges as Ep approaches the Dirac point. 
This is a consequence of the reduced screening of the 
Coulomb interaction. 

It is rewarding to also consider the local TDOS, specifi- 
cally its angular dependence around an impurity. Indeed, 
different kinds of disorder, as parametrized by the ele- 
ments u sl in Eq. (|4]), yield different angular dependences: 
While potential disorder in the Dirac equation yields a 
local TDOS which is rotationally symmetric around an 
isolated impurity, gauge-field disorder, be it abelian or 
nonabelian, leads to a dipole-like angular dependence of 
the local TDOS. Indeed, we find that 



8vl A + B sin 2 



C cos 2 i 



r 3/2 



(k F ry 
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with (j) the angle between the vector r and the direction 
of the bond with modified tunneling amplitude. Here, 
A = u 2 zx + u 2 zy + u 2 zz arises from potential scattering, 
while B = u xz + u 2 yx + u 2 yy and C = u 2 xx + u 2 xy + u 2 yz 
arise from the gauge field induced by the impurity. Thus, 
measurements of the local TDOS can be employed to dis- 
tinguish between the random potential and the random 
gauge field induced by an impurity in the Dirac equation. 

This result is illustrated in Fig. [3l Impurities induce a 
gauge field in the Dirac equation if they modify a hopping 
amplitude. The affected bond then defines the preferred 
direction of the "dipolar" pattern of the local TDOS. 




FIG. 3: Angular dependence of the zero-bias anomaly in the 
coarse-grained TDOS around an impurity inducing a fictitious 
gauge field in the Dirac equation. 

Conclusions — A distinctive feature of the two- 
dimensional electron system in graphene is its exposure 



at the surface, making it directly accessible to local- 
probe techniques including STM. We find that such ex- 
periments can be employed to gain information about 
the character of impurities (potential or gauge-field) by 
exploring the angular dependence of the local TDOS 
around a specific impurity. There is no angular depen- 
dence in the contribution Eq. ([7j) due to disorder scat- 
tering alone. Angular dependence of the local TDOS 
emerges in the zero-bias anomaly due to Coulomb inter- 
actions and impurity scattering, and arises from gauge- 
field disorder only. 
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